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m . Abstract 

o, 

CN I This paper investigates a natural generalization of the K-fi fading channel in which the line-of-sight 

5— ( ' 

(-») , , (LOS) component is subject to shadowing. This fading distribution has a clear physical interpretation, 

■^ ■ good analytical properties and unifies the one-side Gaussian, Rayleigh, Nakagami-m, Ricean, K-/i and 

Ricean shadowed fading distributions. The three basic statistical characterizations, i.e. probability density 

function (PDF), cumulative distribution function (CDF) and moment generating function (MGF), of the 

K-/i shadowed distribution are obtained in closed-form. Then, it is also shown that the sum and maximum 



(yj . distributions of independent but arbitrarily distributed K-fi shadowed variates can be expressed in closed- 

^ ■ form. This set of new statistical results is finally appHed to the performance analysis of several wireless 
communication systems. 

in 

m 

T^lj- I Wireless communications, fading Channels, one-side Gaussian, Rayleigh, Nakagami-m, Ricean, 

-Y-\ . K-fj,, Ricean shadowed. 



Index Terms 



I. Introduction 



o3 . The K-fi fading distribution provides a general multipath model for a line-of-sight (LOS) 

propagation scenario controlled by two shape parameters k and fi. Some classical fading distri- 
butions are included in the K-/i distribution as particular cases, e.g. one-sided Gaussian, Rayleigh, 
Nakagami-m and Ricean. In fact, the fitting of the k-jj, distribution to experimental data is better 
than that achieved by the classical distributions previously mentioned. A detailed description of 
the K-fi fading model and other related models such as 7]-fi and a-yU can be found in [l]-[2]. 
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Shadowing can be introduced in a LOS multipath fading model in two basic ways. The first 
way consists on assuming that the total power, associated to both the dominant components and 
the scattered waves, is subject to random fluctuations. The second way relies on assuming that 
only the dominant components are subject to random fluctuations. The first class of composite 
multipath/shadowing models are named multiplicative shadow fading models, while the second 
class of models are named LOS shadow fading models. 

The Ricean shadowed fading model is a LOS shadow fading model that assumes the Ricean 
distribution for the multipath fading and the Nakagami-m distribution for the shadowing. This 
model fits to the land mobile satellite (LMS) channel experimental data [3], and recently, it 
has been shown that provides an excellent experimental fitting to underwater acoustic commu- 
nications (UAC) channels [4]. In addition, the Ricean shadowed fading distribution has good 
analytical properties and all its basic statistical characterizations are given in closed-form, i.e. 
the probability density function (PDF) and moment generating function (MGF) are given in [5] 
and the cumulative density function (CDF) in [6]. 

Since the k-jjl distribution includes the Ricean distribution as a particular case, a natural 
generalization of the k-ji distribution can be obtained by a LOS shadow fading model with the 
same multipath/shadowing scheme used in the Ricean shadowed model. This paper investigates 
this new model, the k-ji shadowed distribution, so called by analogy with the Ricean shadowed 
distribution. It is shown in this paper that the k-ji shadowed distribution has three main properties: 

• It is motivated by a clear underlying physical model. 

• It provides a remarkable unification of popular fading models: one-side Gaussian, Rayleigh, 
Nakagami-m, Ricean, k-ji and Ricean shadowed. 

• It has good analytical properties; its PDF, CDF and MGF are obtained in closed-form. The 
statistics of the sum and maximum distributions can also be expressed in closed-form. 

The remainder of this paper is organized as follows. In Section II the k-jjl shadowed underlay- 
ing physical model is described. The PDF, CDF and MGF of the k-jjl shadowed distribution are 
derived in Section III, while the sum and maximum distribution are investigated in Section IV. 
Some applications of this set of statistical results to the performance analysis of wireless commu- 
nication systems are presented in Section V. Finally, some conclusions are given in Section VI. 
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II. Physical Model 

The fading model for the k-jj, shadowed distribution relies on a generalization of the physical 
model corresponding to the k-jj, distribution [1]. We consider a signal structured in clusters of 
waves which propagates in a nonhomogeneous environment. Within each cluster the multipath 
waves are assumed to have scattered waves with identical powers and a dominant component with 
certain arbitrary power. While the intracluster scattered waves have random phases and similar 
delay times, the intercluster delay-time spreads are considered relatively large. In contrast with 
the K-fi model which assumes a deterministic dominant component within each cluster, the n-ii 
shadowed model assumes that the dominant components of all the clusters can randomly fluctuate 
as a consequence of shadowing. 

From the physical model for the k-ji shadowed distribution, the signal power W can be 
expressed in terms of the in-phase and quadrature components of the fading signal as follows 

n n 

W = Y,{X, + ip,f + 5^ (r, + eg.)^ (1) 

j=l 1=1 

where n is a natural number, Xi and Yj are mutually independent Gaussian processes with 
E[Xj] = EfYj] = 0, E[Xf] = E[F^^] = cr^, pi and qj are real numbers, and .^ is a Nakagami-m 
random variable with shaping parameter m and E[,^^] = 1. The interpretation of (1) is the 
following. Each multipath cluster is modelled by one term of the sum; thus, n is the number of 
multipath clusters. The scattered components of the ith cluster are represented by the circularly 
symmetric complex Gaussian random variable Xj + jYi. For each cluster the total power of the 
scattered components is 2o-^. The dominant component of the ith cluster is a complex random 
variable given by ^pi + jiqi, thus, its power is given by pf + gf . All the dominant components 
are subject to the same common shadowing fluctuation, represented by the power normalized 
random amplitude ^. 

III. Fundamental Statistics 

This section includes the derivation of the PDF, CDF and MGF of the k-jjl shadowed distribu- 
tion. For the sake of brevity, both the distribution of the signal envelope and the distribution of 
the signal power (or equivalently the instantaneous signal-to-noise ratio (SNR)) will be named 
K-n shadowed; both distributions are connected by a quadratic transformation. 
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The model represented in (1) implies that the conditional probability of the signal power W 
given the shadowing amplitude ^ follows a k-jj, distribution with PDF [1] 



n-l 



1 f w \ '^ _w+d^ ^ I id ^\ 

where (P = J2^=iP'i~^Qi represents the mean power of the dominant components and /j,(-) is the 
modified Bessel function of the first kind [7]. As noted in [1], the natural number of clusters n 
can be replaced in (2) by the nonnegative real extension /i, resulting in a more general and 
flexible distribution. The k parameter is then defined as k = (f/{2a'^fi) and can be interpreted, 
when /i is a natural number, as the ratio between the total power of the dominant components 
and the total power of the scattered waves. In many practical analyses, the random variable 7 
representing the instantaneous SNR is used to model the fading channel; thus, hereinafter we 
will consider the random variable 7 = ^W/W, where 7 = £[7] and W = E[W] = cP + 2cr^/x. 
In terms of the scaled random variable 7, the conditional PDF in (2) can be rewritten as 



Let 7 be a k-ji shadowed random variable with mean 7 and real nonnegative shaping param- 
eters K, fi and m. This fact is expressed symbolically as 7 ~ <Sk^{'J] K,ii,m). The PDF of the 
K-n shadowed distribution is obtained from (3) as follows. 

Lemma 1: Let 7 ~ S^fii'j', k, /x, rn), then, its PDF is given by 

^ r (yu)7(/i/t + m)™ \7/ \ ' ' fj,K + m 7/' 

where i-Fi(-) is the confluent hypergeometric function [7]. 

Proof: See Appendix I. ■ 

It can be checked that the PDF in Lemma 1 is the Ricean shadowed fading PDF when fi = 1. 
Next, the MGF of the k-ji shadowed distribution is derived from its PDF. 

Lemma 2: Let 7 ~ 5k^(7; k, /x, m); then, its CDF is given by 

'J^^ [IJ,K. + m) /„ m(i+«) m ^ 
Proof: See Appendix II. 
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The CDF of the k-jjl shadowed distribution can be expressed in closed-form by the bivariate 
confluent hypergeometric function $2(-)' defined in [7]. 



Lemma 3: Let 7 ~ 5k^(7; k, fi, m); then, its MGF is given by 

/i^-im™ (1 + k)^ (l^^ 



F.(7) = :...::..r,„:i (^) rx 



(6) 

/i(l + /t)7 /i(l + /t) m7 



7 7 fiK, + m^ 

Proof: Taking into account that -^^(7) = C~^[Aiy{—s)/s; 7], the result follows from (27) 

and the Laplace transform pair given in [8, eq. 4.24.3]. ■ 

The fundamental statistics presented in Lemmas 1, 2 and 3 provide a unification of a variety 
of important fading distributions. Table I reflects the parameter specializations which allow us to 
obtain the one-side Gaussian, Rayleigh, Nakagami-m, Ricean, k-jj, and Ricean shadowed fading 
distributions, from the three shaping parameters of the k-ji shadowed distribution. The PDF 
given in Lemma 1 is plotted in fig. 1 for different parameter combinations; it is clearly shown 
the flexibility of the mathematical model represented by the expression (4). 

IV. Sum and Maximum Distributions 

In this Section the distribution of the sum and the maximum of independent non-identically 
distributed (i.n.d) n-fi shadowed random variables are derived. These results indicate that this 
new distribution has good analytical properties, and has great potential as a tool for modelling 
and analyzing a variety of wireless communication systems. 

A. Sum distribution 

The sum distribution of random variables representing the SNRs in a fading channel plays a 
prominent role in the analysis of diversity systems and space-time coding. In the next Proposition, 
the sum of independent k-h shadowed random variables is statistically characterized. 

Proposition 1: Let 7^ ~ S^fii'^k', f^k, /J'k, rnk) for A; = 1, ..., M, where all the random variables 
are arbitrarily distributed and mutually independent. The PDF of the sum 7 = X]a:=i 7fc ^^ given 
by 
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/ 



A (7) 



\ 



1 



V 



r ( E /"fc 1 ^ 

fc=i 



M 

n 



/i^'=mr (1 + n,r [^\ 



^^ ^ -"--^ {fj,kHk + mk)"^'' \-fkJ 



Mfc 



M 



/ 



M 



'±*2 I /^i - "^1, • • • , /^M - "^M, mi, ... , mM] 2_^ ftk] — 



-/iM(l + K,m)1 



k=l 

-fii (1 + Ki) mi7 
7i /ii/ti+mi' 



7i 7m 

-/iM- (1 + k-m) rnMl 



iM 



f^Mt^M + ^M 



(7) 
where $2 (') i^ the confluent muhivariate hypergeometric function [9]. The CDF of 7 is given 



by 



^^7(7) 



/ 



\ 



1 ^ 



/^fc ^k U + /^feJ / J- 



r ( 1 + E /ife 1 ^= 

k=l 



fJ-k 



M 



^(2M) / , , V^ -/il(l + /ti)7 -/iM(l + /tM)7 

$^ M /xi - mi, . . . , /iM - mM, mi, ... , mM; 1 + > Wfc; = , • • • , = , 



-/ii (1 + Ki) mi7 -/iM (1 + k-m) mul 



7i /^i'«i + "^1 



7m /iM/tM + mM 



(8) 



Proof: See Appendix III. 



Once the following technical Lemma is considered, the important independent and identically 
distributed (i.i.d) case for the sum distribution is obtained as a Corollary from the previous 
Proposition. 

Lemma 4: The confluent multivariate hypergeometric function $2 has the following property 



$ 



(AT+M) 



/3i, . . . , /3i, /32, • • • , /32; «^; a;i, . . . , a;i, X2, . . . , X2 = $2 (AA^, (iiM\ v\ xi, X2) , (9) 



N M N M 

where A^ and M are natural numbers, ^[i/] > 0, ^[xi] < and ^[x2] < 0. 

Proof: See Appendix IV. ■ 

Corollary 1: Let 7^ ~ 5k^(7; k, /i,m) for A; = 1,...,M, i.e. all the random variables are 
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^M 



identically distributed and mutually independent. The PDF of the sum 7 = X]a:=i Ik is given by 

1 ^i^^^m^^il + nY^ fiy"^ 



fAl) 



r {M^) (^K + m) 



mM 



] 7^^^-ix 



$2 LuM - mM, mM; /xM; 



— /i(l + /t)7 — /i(l + /t) 7717 



7 



7 fiK + m, 



(10) 



The CDF of 7 is given by 

1 /i^^m"^^(l + /€r'^ /1\^^ 



^^7(7) 



r(l + M/i) (/x/t + m) 



mM 



-1 7^'^x 



$2 /iM - mM, mM; 1 + fiM; 



-/i(l + /T;)7 — /i (1 + /t) 7727 



(11) 



7 7 fiK + m, 

Proof: This result is a direct consequence of Proposition 1 and Lemma 4. 



The asymptotic behavior of the PDF and CDF of the k-(i shadowed distribution is summarized 
in the following result. 

Corollary 2: Let 7^ ~ S^fii'lk] f^k, f^k^'^'k) for k = 1, ..., M, where all the random variables 
are arbitrarily distributed and mutually independent. The asymptotic behavior of the PDF of the 
sum 7 = Yjk=i 7fc when 7;^. ^ 00 for all k is given by 

/ \ 

^ ^ li ifik^k + mkr" V7fc. 



A (7) 



V 



r ( E/ifc 

k=l 



M 



(12) 



/ 



and the asymptotic behavior of the CDF is given by 

/ 

1 



^^7(7) 



\ 



M 



V 



r ( 1 + E /ifc 

fc=i 



n /ifc ^k U + ^fcj /_^\ 
. ^^^ (^,«:, + ,„fc)"^'= V7J 



o/fc — 1 



(13) 



/ 



Proof: This result is a direct consequence of Proposition 1 and the following trivial fact 



AN) 



<l>r'(/3i,...,/3iv;^/;0,...,0) = L 



B. Maximum distribution 

The statistical characterization of the maximum of M independent k-ji shadowed random 
variables is straightforward from the previous results. In general, the CDF and the PDF for such 
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maximum are respectively given by 

M 

''=' (14) 

/max{7,} (7) = ^ n ^-^fe (^) = zi f^k (^) n ^-^^ (^)' 

^ fc=l fc=l r=l,r5>^fc 

where /^^ and F^^ for A; = 1, ..., M are the corresponding marginal PDFs and CDFs. Substitution 
of the expressions for such marginal distributions derived in Section III in (14) provides closed- 
form expressions for the PDF and CDF of the maximum of independent k-ji shadowed random 
variables. 

V. Performance Analysis of Wireless Communication Systems 

This Section shows that the k-ji shadowed distribution is an useful tool for modelling and 
analyzing wireless communication systems. 

In previous Sections it was proved that the k-ji shadowed fading model is a natural general- 
ization of the K-fi model and unifies a variety of popular fading models. Since the k-jjl shadowed 
fading model has an additional parameter m with respect to the k-ji model which is physically 
related to shadowing; the fitting of experimental data to the k-ji shadowed model must be as 
least as good as the fitting to the k-ji model. Otherwise, the same statement is applicable to the 
Ricean shadowed model due to the k-ji shadowed model has an extra shaping parameter ji with 
respect to the Ricean shadowed model. Both the k-ji model and the Ricean shadowed model have 
been proved very useful to model fading scenarios as diverse as mobile radio communications, 
land mobile satellite communications and underwater acoustic communications [l]-[4]; thus, the 
K-fi shadowed model which encompasses these two models represents a very general tool to 
characterize fading channels. 

With regard to the utility of the k-jjl shadowed model for the analysis of wireless communi- 
cation systems, we will show below that the closed-form statistics derived in previous Sections 
allows us to obtain closed- form expressions for certain fundamental performance metrics. In 
particular, the outage probability and/or the error probability for k-ji shadowed fading channels 
will be obtained when the receiver performs maximal ratio combining (MRC) or selection 
combining (SC). These new expressions generalize all the results found in the literature for 
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the K-fi fading distribution and the Ricean shadowed distribution, and all the fading distributions 
encompassed by these two models. 

A. Selection combining with k-ji shadowed fading 

Let us consider a receiver with L branches which performs SC. Each branch experiences 
K-fj, shadowed fading with an instantaneous SNR 7^ ~ ^K^i'lk', f^k, f^kj'm'k) for k = l,...,L. It 
is assumed that all the random variables 7^ are mutually independent. Then, using (14) and 
Lemma 3, the outage probability for SC is given by 

Mfc 



^^^ r (/ifc) (/ifc/tfc + mfc) V7fc 



ri^'^x 



(15) 

^ , /ife(l + /tfc)r7 ;Ufe(l + /tfe) murj 



Ik Ik fikKk + rrik, 

where 1] is the SNR threshold. Since the function $2 tends to unity when 7fe — )> for all 
k = 1,...,L. After taking into account that $2 (/^i^ •••^/^at; '^; 0, ..., 0) = 1, the asymptotic 
behavior of Pq is given by 

fJ-^T {nk) ifikKk + mk) X'jkJ 
Fig. 2 shows the outage probability for SC computed by (15), and superimposed simulation 
results which validate the analytical derivations. Some comments on the numerical computation 
of the $2 function are presented in Appendix V. In Fig. 2 it is assumed a particular scenario 
with three branches for SC in which 71 = 72 = 73 = 7, ^i = 1-2, K2 = 2.7, k^ = 3.1, /ii = 4, 
/X2 = 2, /i3 = 1 and mi = m2 = m^ = m. The curves represent the outage probability in terms 
of the average SNR per branch 7 for different values of the shaping parameter m. The results 
for this particular scenario show the significant impact of shadowing in the system performance, 
despite the k parameter which measures the LOS strength is below 5 dB at every branch. When 
m — 7- 00 these results are showing the performance of SC when fading is of k-jjl type. 

B. Maximal ratio combining with n-fi shadowed fading 

In this subsection we consider a receiver with L branches which performs MRC. Each branch 
experiences n-fi shadowed fading with an instantaneous SNR 7^ ~ 5K^(7fc; Kfc, Atfc,^A:) for 
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k = I, ..., L.ltis assumed that all the random variables 7^ are mutually independent. The outage 
probability is straightforward from Proposition 1 

/ \ 

1 



Po = Pr {-fMRC ^ V} 



V 






k=l 






$ 



(2L) 



Hi - mi, . . . , fXL - rriL, mi, . . . , rriL 



i + Z 



^fc; 



-/ii(l + fi;i)?7 -/ii(l + Ki)?7 



fc=i 
-;Ul (1 + Ki) miT] 

7i /iiKi + mi' 



71 



-fJL (1 + /t^) 

7l 



7l 

17111] 



fXLK,L + rriL 



(17) 



where 77 is the SNR threshold. The asymptotic behavior of the outage probability when 7fc — > 
for all k = I, ..., L is directly obtained from Corollary 2 

/ \ 



Pn 



1 ^ fit'K' (1 + '^kf n^^' 



r ( 1 + E /"fc 

k=l 



^ ^ fcJi il^kK,k + rrikf^'' \-ik 



J2 f^k 
/rtk—l 



(18) 



/ 



Now we will prove that the bit error probability of MRC systems under K-fi fading can be 
computed in closed-form. The bit error probability of many wireless communication systems 
with coherent detection is determined by 



R 



n = 5^«rE Q UA 



r=l 



r7 



(19) 



where {a^, (3r}^=i are modulation dependent constants [10]. For MRC, the bit error probability 
can be obtained from (20) after integrating by parts. 

R 



-"■ poo 

Pb = ^ar Q (^v^j fiMRc (7) d-f 

R 

E 



O^rVPr I 6 2' 



(20) 



IMRC 



(7) dl- 



,=1 vSvr io v^ 

Substituting (6) in (20) and using [9, pp. 290, eq. 55], the following closed-form expression is 
obtained 
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\ 



L 



8^ \J. 



r=l 



X 



L L 

E "^ 



/ifc, fii - mi, . . . , iJ,L - rriL, mi, ... , m^,; 1 + 2^ Atfc; 



-2/ii (1 + Ki) 



(21) 



-2/ii (1 + «l) -2/Ui (1 + Ki) mi -2/ii (1 + /^l) m^ 



lL(3r llPr /Ul/tl + mi ^iPr I^Lf^L + ^L 

where Fj^ (■) is the multivariate Lauricella function [9]. 

The outage probability for MRC computed by (17) is plotted in fig. 3, including superimposed 
simulation results which validate the analytical derivations. The numerical computation of the 
multivariate $2 function is discussed in Appendix V. The same particular scenario used for 
fig. 2 is assumed here; i.e. MRC with three branches in which 71 = 72 = 73 = 7, f^i = 1-2, 
K2 = 2.7, /t3 = 3.1, /ii = 4, /i2 = 2, /i3 = 1 and mi = m2 = m^, = m. In this figure, the outage 
probability for MRC is plotted as a function of the average SNR per branch 7 for different 
values of m. As in the SC case, the shadowing parameter m has a great influence on the system 
performance. 

The bit error probability for MRC is plotted in fig. 4 when a BPSK modulation is used, 
i.e. i? = 1, tti = 1 and (3i = 2. Fig. 4 displays both analytical results computed by (21) and 
simulation results. The numerical computation of the multivariate F^ ' function is discussed in 
Appendix V. Again, the same particular scenario used for fig. 2 and fig. 3 is assumed here. In 
this figure, the bit error probability for BPSK with MRC is plotted as a function of the average 
SNR per branch 7 for different values of m. As with the outage probability, the shadowing 
parameter m has a great impact on bit error probability. 

VI. Conclusions 

The statistics of the n-fi shadowed fading model have been derived along this paper. This 
fading distribution is a natural generalization of the k-jj, fading channel which includes shadow- 
ing. Such fading distribution has a clear physical interpretation, good analytical properties and 
unifies the one-side Gaussian, Rayleigh, Nakagami-m, Ricean, /t-/i and Ricean shadowed fading 
distributions. The three basic statistical characterizations, i.e. probability density function (PDF), 
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cumulative distribution function (CDF) and moment generating function (MGF), of the /t-/i 
shadowed distribution are obtained in closed-form. It is also shown that the sum and maximum 
distributions of independent but arbitrarily distributed k-jj, shadowed variates can be expressed in 
closed-form. The derived closed-form statistics are given in terms of the bivariate hypergeometric 
functions $2 and Fu or the multivariate functions $2 ^"^^ ^d • Numerical methods to compute 
these functions have been discussed. Finally, this set of new statistical results is applied to 
the performance analysis of several wireless communication systems. In particular, the outage 
probability and the bit error probability for systems employing SC and MRC over k-ji shadowed 
fading channels have been investigated. 



Appendix I 
Proof of Lemma I 



From (3), the PDF of 7 can be computed as 



00 



A (7)=/ kdr,i)fdi)di 
'0 



' > e 7 — - — -0 (7] 



M-i 



(22) 



7/t^ v7/ r("^ 

where 



r 2e-«^(M-+™)^2„.-^j^_^ I 2/ieW ^^^^^ I di. (23) 



9(7) 

Jo 

The quadratic transformation (t = ,^^) in the integral which appears in 0(7) yields 



6(7) = yJ|°'e-(-+-)r-t-^V, {2^^J^^±^ 



t dt. (24) 



Sequential application of the identities [8, eq. 4.16.20] and [7, eq. 9.220-2] allows us to express 
6(7) in terms of the confluent hypergeometric function iFi 

9(7) = rM(^^^(l + ^)r^ (A^.F, L,/.;^^Mi±^2^ . (25) 

r (/i) (nK + m)"' \7/ \ ' ' fiK + m 7/ 

The proof is completed after plugging (25) in (22) and performing some algebraic simplifications. 
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Appendix II 
Proof of Lemma II 

Taking into account the linearity and the frequency shifting properties of the Laplace transform 
yields 



M,{s) = C[f,{^);-s] 



fi'^m"' (1 + Kf 
r (/i) {fXK + m)"* ^7 



X 



c 



Y-\Fi(m,fx 



/i^/t(l + k) -f\ ji {1 + k) 
IXK + m 7/ ' 7 



— s 



(26) 



The Laplace transform in (26) is recorded in [8, eq. 4.23.1]; thus, the MGF can be expressed as 



M {s)= -- .^...^ (lyrM A_/i(l+i£ 

''' r (/i) (/i/t + m)™ \7/ (-s)^ V -57 



/i (1 + k) m 
S7 fiK, + m 



(27) 



which after some straightforward algebraic manipulations takes the form expressed in (6). 



Appendix III 
Proof of Proposition I 



The MGF of the sum distribution is given by 



M 



M,{s) = l[ 



f^k ^k U + '^k) / i 



\nife 



\ {l^kK-k + mk) - V7fe 



Mfe 



1 / , _ /ifc (1 + K,k 



From (28), the PDF of the sum can be expressed as 

/ 



f,{^) = C~'[M,{-s);^] 



12k (1 + Kk) rrik 
sik l^kK.k + rrik 

\ 



-nik 



V 



M 



\Mfc 



M 

r ( E /"^ "== 

A:=l 



;^ (/ifc/tfc + "^fe)""° V7fc 



Mfe 



X 



/ 



c 



-1 






(28) 



M 

E Mfe fc=l 



n 1 



k=l 



12k (1 + /tfc) rrik 
sik f^kK,k + rrik 



-nik 



7 



(29) 



April 30, 2013 



DRAFT 



14 



In such arrangement, the right side of (29) can be identified with [9, pp. 290, eq. 55] yielding 
the expression for the PDF stated in the proposition. To obtain the CDF we can observe again 
that 

/ \ 



F,{^) = C~^[M,{-s)/s-^] 



1 ^ 



^l^^''wJ^'' (1 + KkY" ( 1 



r ( 1 + E /"fc 1 '^^ 

fc=i 



\ ifj'k'^k + mfc)™*= \jk 



t^k 



X 



/ 



C 



-1 



M 



k=l 



/ifc (1 + Kfc) ruk 



S * = ! 



-nife 



i7 



(30) 



A new identification of (30) with [9, pp. 290, eq. 55] completes the proof. 



Appendix IV 
Proof of Lemma IV 



Let us consider the following ancillary function 



A (t) = t'^-l.l.f +^) /32, . . . , /32, /32, . • • , /32; l^] Xit, ..., X,t, X2t, ...,X2t 



(31) 



N 



M 



N 



M 



Considering the Laplace transform of A (t) which is obtained with the help of [9, pp. 290, eq. 55], 
performing trivial algebraic simplifications in the transformed domain, and returning again to 
the t-domain with [9, pp. 290, eq. 55] yields the required property after setting t = 1. 



Appendix V 
Numerical computation of the functions $2 and F^ 

Most of the results derived in this paper involve either the bivariate functions $2 and Fjj or the 
multivariate functions $2 ^^'^ ^d ■ Therefore, some comments on the numerical computation 
of these special functions can be useful for the reader. Each of them will be treated separately. 
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The bivariate hypergeometric function Fd is the same as the Apell hypergeometric function 
Fi, which is implemented in the most popular scientific software packages, e.g. MATLAB and 
MATHEMATICA. Therefore, its computation is straightforward by these software tools. 

The bivariate confluent hypergeometric function $2 is defined in the popular mathematical 
handbook edited by Gradshteyn and Ryzhik; however, it is not yet implemented in MATLAB and 
MATHEMATICA. As with the Marcum Q function which has a Bessel series representation, the 
$2 function can be expressed as a iFi series which is very appropriate for numerical computation 
[11, eq. 4.19] 

$2(^ b'; c; w;z) = f] ^^w\F,{h'- c + k;z). (32) 

The multivariate hypergeometric function F^ is not yet implemented in MATLAB and 
MATHEMATICA; however, it can be easily computed by its Euler-type representation and 
standard numerical integration methods 

Fjf)(a,6i,...,6^;c;a;i,...,a;;v) = ttt^ V / t^-\l-tr"'~\l-x^tr''\..il-XJ,t)-''^dt, 

r(a)r(c-a) Jq 

(33) 
where ^[c] > ^[a] > 0. Note that this last condition is satisfied in the multivariate F^ ' function 
which appears in (21). 

The multivariate confluent hypergeometric function $2 is not yet implemented in MATLAB 
and MATHEMATICA; however it can be efficiently computed by inverting its one-dimensional 
Laplace transform [9, pp. 290, eq. 55]. Numerical methods for inverting Laplace transforms are 
exhaustively discussed in [12]. 
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TABLE I 

Common Fading Distributions Derived from the k-^j, Shadowed Distribution 



Fading Distribution 




Parameters of the k-^ Shadowed Distribution 


One-sided Gaussian 


fl = 0.5, K ^ 0, TO ^ oo 


Rayleigh 


^=1, K^O, TO^OO 


Nakagami-m, 
with shaping parameter m 


/i — ?7i,, K— >0, ?Tt— >00 


Rician, 
with shaping parameter K 


/i— 1. K — _y. ?7i,— >00 


with shaping parameters k and 


M 


/i = /i, K = K, TO— >00 


Rician shadowed, 
with shaping parameters K and 


m 


fi= 1, K = K, m = m 
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PDF 



2r 



Rician Shadowed 
K=5,/i=l,m=2 




Fig. 1. PDF of the k-^ shadowed distribution (7 = 1). 
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Fig. 2. Outage probability versus average SNR per branch over n-fj, shadowed fading channels. A triple-branch SC scenario 
is considered, with parameters 71 = 72 = 73 = 7, ki = 1-2, H2 = 2.7, K3 = 3.1, fii = 4, ^2 = 2, fj, = 1 and 

mi = 7712 = mi = 771. 
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Fig. 3. Outage probability versus average SNR per branch in K-fi shadowed fading channels. A triple-branch MRC scenario 
is considered, with parameters 71 = 72 = 73 = 7, ki = 1-2, K2 = 2.7, res = 3.1, fii = 4, 112 ~ 2, fj, = 1 and 
mi — m2 ~ m-i = m. 
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Fig. 4. Bit error rate versus average SNR per branch in K-fi shadowed fading channels. In this plot BPSK modulation and 
triple-branch MRC are considered, with parameters 71 = 72 = 73 = 7, «;i = 1-2, K2 = 2.7, K3 = 3.1, fii = 4, fi2 = 2, jj, = 1 
and mi = m2 = ma — m. 
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